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The baby Skyrme model is studied with a novel choice of potential, V = This "easy 

plane" potential vanishes at the equator of the target two-sphere. Hence, in contrast to previously 
studied cases, the boundary value of the field breaks the residual SO (2) internal symmetry of the 
model. Consequently, even the unit charge Skyrmion has only discrete symmetry and consists of a 
bound state of two half lumps. A model of long-range inter-Skyrmion forces is developed wherein 
a unit Skyrmion is pictured as a single scalar dipole inducing a massless scalar field tangential to 
the vacuum manifold. This model has the interesting feature that the two-Skyrmion interaction 
energy depends only on the average orientation of the dipoles relative to the line joining them. 
Its qualitative predictions are confirmed by numerical simulations. Global energy minimizers of 
charges B = 1, . . . , 14, 18, 32 are found numerically. Up to charge B = 6, the minimizers have 
2B half lumps positioned at the vertices of a regular 2B-gon. For charges B > 7, rectangular or 
distorted rectangular arrays of 2B half lumps are preferred, as close to square as possible. 

PACS numbers: 11.27.+d, 05.45.Yv, ll.10.Lm, ll.10.Jj 



I. INTRODUCTION 



The Skyrme model, and its smaller relative, the baby 
Skyrme model, have received much attention during the 
last few decades. The models are very interesting from 
the mathematical physics point of view: they contain 
topologically nontrivial field configurations called topo- 
logical solitons and the Skyrme model arises as a large- 
TV limit of QCD 0, @|. The baby Skyrme model is the 
two-dimensional analog of the full model and its main 
interests lie in its easier tractability, the possibility of 
using the baby Skyrme model results as guidelines in in- 
vestigating the full Skyrme model and some condensed 
matter contexts, where the baby Skyrme model itself is 
of physical relevance Q. The baby Skyrme model was 
investigated intensively during the 1990s but interest in 
the model has recently received a new boost [3-Hlj. ap- 
parently due to the prevalence of baby Skyrmions in var- 
ious condensed matter systems [12, sucn as Fei-^Co^Si 
[l3|, MnSi flil . [l5j and quantum Hall systems (for a re- 
view, see [lj| and references therein), despite the fact 
that these are not modeled by the baby Skyrme model. 
It is against this background that we investigate the baby 
Skyrme model. 

The behavior of the solitons of the baby Skyrme model 
depends crucially on the choice of potential. For some 
potentials, there will be no smooth solutions at all; for 
some, like that studied in [l?}, there will be exact solu- 
tions for the equations of motion and for the rest, such as 
here, numerical methods must be used. So far, only cases 
where the potential and boundary condition preserve the 
SO (2) symmetry of the Lagrangian have been studied in 
detail. Here we study the easy plane potential, which 



breaks 0(3) symmetry to SO(2), and the boundary con- 
dition further breaks the symmetry to This endows 
the model with one massive and one massless mode and 
one expects significantly different behavior from 5*0(2) 
symmetric cases. 

We use numerical methods to demonstrate that the 
model has a rich static energy structure, including lo- 
cal minima and that multi- Skyrmions are stable in this 
model. We also develop an approximate way to estimate 
the long-range interaction of two unit baby Skyrmions 
and verify this numerically. We identify some open ques- 
tions and present a conjecture about the energy of the 
global minimum, when the charge approaches infinity. 



II. THE MODEL 



The Lagrangian density of the model, defined on 



is 



C = \c 2 d^ a d^ a - \c A F^F^ - CoW^) 2 , (1) 



F u 



(2) 



where the first terms are called the Dirichlet, the Skyrme 
and the potential terms, respectively and Cj are free pa- 
rameters (coupling constants) of the model. The sub- 
scripts are chosen to match the number of derivatives in 
the term and will prove useful later. 

Choosing the usual metric (+,—,—) yields the static 
energy density 



£ =|c 2 ||V« 



(3) 



If one requires the solitons to have finite energy, it is 
necessary to impose the boundary condition 
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and choose from the vacuum manifold of the model, 
i.e. such that 

17(6x0 = 0. (5) 

As usual, Eq. (0| allows us to consider the field as map 
: S 2 — > S 2 , which gives the system its interesting topo- 
logical properties, namely, its fields arc classified topo- 
logically by the homotopy group 7r 2 (S 2 ). The homotopy 
invariant of the field is called the (topological) degree and 
is given by the expression 

B=-j- [ f (WxW)d 2 i. (6) 
4?r J R 2 

Furthermore, one can expect the system to support 
topological solitons based on two factors. First, Derrick's 
theorem of nonexistence does not apply since the terms 
zcroth and fourth order in derivatives behave in oppo- 
site ways upon uniform scaling (the second order term is 
scale invariant), thus allowing solutions with a preferred 
size. The second reason is the existence of a Bogomol'nyi 
bound for the (static) energy which depends on the de- 
gree of the field, given by [ill ] 

E > 47r(c 2 + y/wi(U)) \B\ , (7) 

where (U) is the average value of U : S 2 — > M. Unlike in 
some models, the bound cannot generically be saturated 
in this model. 

In this work, we take £/(</>) = (f>3 and (poo = (1,0,0). 
This potential has the exceptional effect of allowing the 
breaking of the symmetry of the energy almost com- 
pletely. While the energy density exhibits SO (2) symme- 
try (from the potential) , the choice of boundary condition 
breaks this down to Z 2 . Hence, even the Skyrmion with 
unit topological charge will have only discrete symmetry. 

We seek static solutions in various homotopy classes 
and investigate the forces between some of these. In order 
to find a static solution, we numerically minimize the 
static energy of an initial field configuration. Some initial 
configurations will be constructed from existing minima, 
but those that are obtained from an exact expression are 
taken to be 

w = X(x + iy) B (8) 

i /M 2 -i\ 

4> = —2— 2ReH , (9) 
M + 1 \-2Inu»/ 

which for B — 1 gives 

where A is simply a scale factor to adjust the size of the 
initial configuration. We also added random ripples to 
destroy any discrete symmetries left over from the ini- 
tial configuration. It was observed that this is usually 



unnecessary, but some minimizcrs take different orienta- 
tions (but keep their shape) when the random ripples are 
omitted. 

In all that follows, we set Co = c 2 = C4 = 1.0. 



III. THE UNIT SKYRMION 

The lattice approximation of derivatives is achieved by 
using a simple forward differencing scheme. The discrete 
energy functional is then minimized using gradient based 
methods. 

The gradients are computed from the discrete total en- 
ergy in a straightforward manner. The energy minimiza- 
tion is then achieved with the TAOQJ] and PETScdHHl 
parallel numerical libraries. The libraries provide sev- 
eral different minimization schemes, of which we have 
chosen the (in this case) fastest algorithm, the limited 
memory quasi-Newton algorithm (also called a variable 
metric algorithm) with BFGS |22l - l25l ] formula for Hessian 
approximations. We consider the minimization to have 
converged, when the sup-norm of the gradient is less than 
10 -7 . This may sound like a very strong requirement, but 
it turns out to be necessary. Our initial choice of 10 
works for most cases, but when the results were checked 
with a the stronger requirement, they experienced sig- 
nificant changes and decreases in energy after attaining 
the looser condition. Using a still stronger criterion did 
not change the energies or shapes of the minimizers any 
further. 

The accuracy of the numerical scheme was tested 
against the known exact solution of w = (x + iy) (26j 
for the potential U = (1 — (fe) 2 . We find that the accu- 
racy of the numerical scheme is heavily dependent on the 
size of the lattice, due to the energy density of the Dirich- 
let term being slowly decaying. When in a large enough 
lattice, the difference in energies between our minimizer 
and the exact solution is less than 1%, so we expect sim- 
ilar accuracy in the present model as well, since it has 
similar decay properties. 

We always use a square lattice and will refer to the 
number of lattice points along its side by a number N . 
Starting with (fTTJj) on a N = 801 lattice with lattice con- 
stant h = 0.1, the corresponding minimum energy den- 
sity is depicted in Fig. Hal The broken rotational sym- 
metry is immediately evident in the form of the two red 
peaks. A view of how the field cf>3 itself behaves is shown 
in[Tb]and an immediate observation is that the half-plane 
y > gets mapped to the upper hemisphere of the target 
S 2 and the half-plane y < gets mapped to the lower 
hemisphere. This, combined with the fact that each in- 
dividual peak in energy density strongly resembles a CP 1 
model lump, leads us to call the solution on a half-plane 
a half lump. 

Note that the energy density image is zoomed in on 
the central area of the numerical lattice, where almost all 
the energy density is concentrated. Obviously, the whole 
computation could have been done in a much smaller 
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lattice, but we opted not to find the optimal lattice sizes 
for each and every configuration due to the fact that the 
amount of computer resources saved would have been 
small. We also confirmed that the resulting field con- 
figuration is not dependent on either the value of the 
lattice constant or "physical" size of the computational 
box by first repeating the computation with lattice con- 
stant h = 0.1 and N = 1599 and then repeating again 
with h = 0.05 and N = 1599. 

It is instructive to look at the fields themselves, in par- 
ticular, their asymptotics at large r. Since — > (1, 0, 0) 
as r — > oo, we must have 



</> = (y/i- <h, <fe) = (i, o, 0) + (0, 2 , 3 ) + ■ ■ ■ 

at large r, where 02,03 are small. Substituting this into 
C and keeping only leading order terms yields the La- 
grangian density of an uncoupled pair of real scalar fields, 



C 



C2 \ 7^02^02 



A A C2 



(12) 



one massless (0 2 ) and the other of mass y/co/c 2 (03). 
One expects the large r field of a B = 1 soliton to be well 
approximated by a static solution of this linear model. 
Since B = 1, 02 + i03 should wind once around as one 
traverses a circle at large r. This leads one to predict 
that 



<12 



<h(r,e) ~ ^costf, 03 (r,6) 
r 



q3Ki{\/c /c 2 r) sin 6» 
(13) 

at large r, where q 2 , q% are unknown constants. 

Numerical evidence in favor of the conjectured asymp- 
totics for is presented in Fig. [51 which shows plots of (gj) 
r02 (r, 9) and jb| 03 (r, 9) /K\ (r) against 9 for an increas- 
ing sequence of values of r. By fitting sinusoidal curves, 
we find that q 2 ~ 3.1 and q 3 rj 9.9. Note that, at large 
r, K\{r) decays like e~ r /\/r, so the massive field 03 de- 
cays much more quickly than the massless field 02. This 
has important consequences for the long-range forces be- 
tween Skyrmions, as we shall see in the next section. 




(a) The energy density of the (b) The contour plot of 03. 
central region of the lattice. 

Figure 1: The minimizcr for degree B = 1. 



Having found the B — 1 solution, it is natural to ask 
whether there are multisolitons, that is localized solu- 
tions with B > 1. Given the linear dependence of the 




(a) Plot of rcf>2(r, 8) against 
and best fit for 92 cos 8. 



(b) Plot of </> 3 (r,fl)/A"i(r) 
against 8 and best fit for 
53 sin 8. 



Figure 2: Plots of 2 and 3 for r e {2, 3, 4, 5, 7.5, 10} 
(the dashed lines, colored blue, green, red, cyan, 
magenta, and yellow, respectively) and the 
corresponding best fits (solid black). 



energy bound on B, Eq. ©, it is conceivable that all 
these would consist of B copies of the B = 1 solution at 
infinite separation. Indeed, this is what happens for the 
potential U = (1 — 03) 2 , where the localized B > 1 solu- 
tions are unstable against decay to B infinitely separated 
solutions with B = 1 [l7| . 

In that model (the so-called "holomorphic" model) , the 
force between a pair of unit charge skyrmions turns out 
to be repulsive, whatever their separation and relative 
orientation. By contrast, we will see that in the current 
model (as for several other choices of potential) there is 
a choice of relative orientation for which solitons attract 
one another at long range (the so-called attractive chan- 
nel), making it plausible that spatially localized solutions 
of every charge exist. 



IV. INTERSKYRMION FORCES 



Intersoliton forces in the baby Skyrme model with an 
easy-axis potential ^[/(0) 2 = 1 — 03 were analyzed in 
detail by Schroers et al in [27}, and their methods read- 
ily adapt to the current setting. We shall assume that 
the dominant interaction at long range is mediated by 
the massless field 02 and ignore the contribution of the 
faster decaying field 03. The key observation is that the 
asymptotic field 



Mr, 6) 



q-i 



cos 9 



(14) 



coincides with the static solution of the linearized model 



£ = ^026^02 + «02 



(15) 



in the presence of an appropriate point source k at the 
origin, in this case 



2tt 



d x 5{x) 



(16) 



where 5 denotes the two-dimensional Dirac delta distri- 
bution. To see this, note that the elementary kernel [28| 
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for Laplace's equation on M 2 is G(x) = (2%) 1 log | a; | , 
that is, 



G, 



S(x) 



(17) 



Source (| 10[) induces the field in standard orientation and 
position. The field rotated (spatially) through angle a 
and translated to position y £ R 2 is induced by the 



k(x) = d ■ WS(x — y) 



(18) 



where d = — (72 (cos a, sin a). Physically, one interprets 
this as the masslcss scalar field induced by a dipole source 
of moment d at position y. It is important to realize that 
the dipole moment which replicates the soliton's asymp- 
totic field points orthogonally to the line connecting its 
constituent half lumps. One should not, therefore, think 
of the half lumps as being scalar monopoles of opposite 
charge, generating a dipole moment (recall that the half 
lumps are structures in the 03 field, while the long-range 
field is dominated by 02). 

Now consider the situation where two B = 1 solitons 
with orientations a±, a 2 are placed at positions y, z with 
R = \y — z\ large. The interaction energy that they 
experience should, on physical grounds, coincide asymp- 
totically with the interaction energy of the corresponding 
scalar dipolcs in the linear theory, that is, 



BU = - ^(x)^\x)d 2 x 



(19) 



where «W(x) = d\-V 5{x—y) and is the field induced 
by dipole d2 at z, that is 



4 2 \x) = - — d 2 -V(log\x-y\). (20) 



Hence 



E int = — [ d 1 -V[S(x-y)]d2- V(log \x - z\) d 2 x 



2tt 



di ■ V y d 2 ■ V z / 5(x - y) log \x - z\ d 2 x 



2nR 2 



{di ■ d-2 — 2(n • di)(n ■ d 2 )} 



(21) 



where n = (y — z)/\y — z\, the unit vector directed from 
dipole 2 to dipole 1. It is interesting to note that this is 
precisely minus the interaction energy of a pair of electric 
dipoles in two-dimensional electromagnctism, which can 
be understood by noting that, in scalar field theory, like 
monopole charges attract, while opposite charges repel. 
In the case where the solitons are placed symmetrically 
about the origin on the x\ axis, that is, y = (R/2,0), 
z = (—R/2,0), one sees that 



9 

2ttR 2 



cos(ai + a 2 ) 



(22) 



This formula is rather remarkable. It has the interesting 
property that it depends only on the average orientation 



(a) 



attract 



(c) 



repel 



■1.4 



attract 



(d) 



Figure 3: Predicted behavior of two-Skyrmion 
superpositions. The disks represent half lumps with 
03 > (red) and 3 < (blue), while the black arrows 
represent the scalar dipole moment associated with each 
Skyrmion. 



Oil + et 2 of the dipoles, and is independent of their phase 
difference a\ — a 2 . This is precisely the opposite of the 
situation found by Schroers et al in the easy-axis model 



E, 



4 



easy — axis 



— K\(R) cos(ai 
2tt 



"2), 



(23) 



which highlights a fundamental difference between the 
two models. The unit soliton in the easy-axis model is 
axially symmetric, and spatially rotating the soliton is 
identical to internally rotating it (about the axis 0oo)- 
Consequently, the interaction energy in the easy-axis case 
must be invariant under (ai,a 2 ) H> [a\ + A,a 2 + A), 
and hence can depend only on oc\ — a 2 . In the easy plane 
model, however, the unit soliton has no rotational sym- 
metry. One cannot rotate it internally about 0^ (the 
fields so obtained do not solve the field equation), and 
spatially rotating the soliton produces a genuinely phys- 
ically distinct solution (e.g. its energy density changes). 
So there is no reason to expect Ei nt to be independent of 
0:1 + That it depends only on a\ + a 2 is more sur- 
prising, and is a basic fact about dipole interactions in 
two dimensions. In practice, it seems likely that higher- 
order corrections will break the independence of Ei nt on 
a.\ — a 2 . 

The most important consequence of formula (f2"2"| is 
that it predicts long-range attraction between unit soli- 
tons if they are appropriately oriented relative to one an- 
other, leading one to expect there to exist stable, static, 
spatially localized solutions in all homotopy classes. In 
particular, soliton pairs initially oriented as in Fig. [3ta) 
oi'EHtj) are predicted to attract one another, while those 
oriented as in Fig. (3Jc) or [2Jd) are predicted to repel. 
These predictions were vindicated by the numerical sim- 
ulations described in the next section. 
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V. MULTISKYRMIONS 

In this section, we report the numerical results on 
multi-Skyrmions, confirming the validity of the dipole ap- 
proximation, and give examples of the very complicated 
landscape of local minima of higher multi-Skyrmions. 
We split this section in two according to two different 
"branches" in the shapes of the minima. Before going 
further, it is worth noting several general observations. 

First, the dipole approximation was confirmed for pairs 
of B = 1 solitons. To confirm this, we used three meth- 
ods: the usual BFGS method, simple gradient flow and 
also a Newton flow second order dynamics, where we set 
ip := do4> and do4' = —d^E, do4> = ip. They all agree that 
two unit Skyrmions starting in the attractive channel co- 
alesce into a B = 2 multisoliton with an energy lower 
than twice that of the unit Skyrmion. The resulting two- 
Skyrmion is depicted in Fig. | ^a)| However, confirming 
the repulsive channel is slightly more challenging. First, 
the computational box needs to be large enough, other- 
wise the boundaries will exert pressure on the solitons 
and they react by rotating into the attractive channel 
and coalescing. Second, the solitons need to be placed 
exactly symmetrically on the lattice, lest they again ro- 
tate into the attractive channel and coalesce. When set 
up correctly, however, placing two B = 1 Skyrmions in 
the repulsive channel results in them receding. Thus we 
can confirm the qualitative correctness of the dipole ap- 
proximation. 

Second, all the results presented below are done in a 
standard lattice with ft. = 0.1 and TV = 801, except the 
case B = 32, which docs not fit into the smaller lattice at 
all. For those values of B, where multiple minima were 
found, two additional methods were used for their ini- 
tial configurations. The first method is a simple cut and 
paste of other solitons, for example, two 2-solitons were 
put side by side to produce a 4-soliton. The other method 
was to start with TV = 801, h = 0.05 which resulted in a 
different minimum than TV = 801, h = 0.1. This is due 
to boundary pressure: if TV was doubled, the other result 
reemerged. When this happened, the h = 0.05 solution 
was made sparser by a factor of 4 to put it into the stan- 
dard lattice, where it was then minimized again to get 
a result which can be compared with the others. Where 
several initial states produced the same minimizer, we 
only mention one of the initial states. 

A. Charges B £ [2,6] 

Minimal energy Skyrmions of charges B in the inter- 
val [2,6] found using our numerical scheme form regular 
polygons, where there are half lumps at the 2B vertices. 
The energies of these minimizers are always lower than 
BEi, where E\ is the energy of the unit Skyrmion, so 
there is good reason to believe they are stable in ac- 
cordance with the above dipole picture. It should be 
noted, however, that the symmetry of these minimizers 



is the same as that of the initial data w = X(x + iy) B . 
Therefore, one is led to suspect that the random ripples 
added to w may not be enough to break the symmetry. It 
is therefore necessary to introduce a stronger symmetry 
break. For B = 2, the computations done to confirm the 
dipole approximation provide this break, and still relax 
to the same square final configuration. The obvious way 
to test the higher charges was therefore to place two lower 
charge solitons side by side and minimize their energy. 

Indeed, for charges 5 and 6, there are minima, which 
were only found using this spliced initial configuration. 
However, they have higher energy than the respective 
polygons, and therefore represent a local minimum only. 
For the charges 3 and 4, various different initial configu- 
rations were tried, but they always relaxed to the respec- 
tive polygon. We present the minimizer of B = 6 as an 
example of the polygon shaped minimizers in Fig. | ^b)| 

B. Charges B > 6 

The 273-gon form of the global minimizers only persists 
up to B = 6 after which, the shape of the minimizer is 
completely different and there is no obvious rule for the 
shape. We will briefly discuss these minimizers below, 
but first it should be noted that the symmetric (plus ran- 
dom ripples) initial state still produces a polygon shape 
local minimizer for all charges except 13 and 14, but this 
is no longer a global minimum. The polygon exists for 
charges 13 and 14, also, but was only obtained by start- 
ing the minimization process with TV = 399, h = 0.1 and 
embedding the result in the standard lattice. The poly- 
gons are, however, quite stable and the energy differences 
between them and the global minima are very small. For 
example, we squeezed some of the polygons to half the 
size in the y direction and they relaxed back to the poly- 
gon shape instead of the global minimum. 

The solutions obtained from lower charge solitons set 
side by side and starting the minimization in a TV = 
801, h = 0.05 lattice, produce a veritable zoo of different 
minimizers, of which all but one have lower energies than 
the 273-gons, the exception being charge 14. Because 
of the amount of work involved in patching the solitons 
together and the total number of computations already 
performed, we did not attempt an exhaustive search with 
side-by-side-type initial configurations of all charges. 

The first nonpolygon minimizer occurs at B = 7 and 
looks like a sligh tly d istorted rectangle of 2 x 7 half-lumps. 
It is depicted in^c)]and was obtained by putting a B = 3 
soliton next to a B = 4 one. A similar minimizer is 
also found for B = 9 using a B = 4 soliton beside a 
B = 5 soliton as the initial configuration. These are 
global minima, within the limits of our numerics. 

At B = 8 the half lumps of the minimum energy config- 
uration arrange themselves into a square lattice of 4 x 4 
half lumps, as seen in Fig. | ^d)[ The initial configura- 
tion consisted of four B = 2 solitons at the vertices of 
a square. Similar square lattice configurations are also 
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found for B = 18 and B = 32, shown in Figs. $ej\ 
and[ ^f)| The initial configurations were as follows. For 
B = 18, two B — 4 and two 5 = 5 solitons were put at 
the opposite vertices of a square and for B = 32, four 
B = 8 solitons were placed at the vertices of a square. It 
is interesting to note that the B = 2 solution can also be 
thought of as 2 x 2 square lattice of half lumps, providing 
us with a tantalizing opportunity to conjecture that all 
2B = (2q) 2 ,q € Z will fall into similar square lattices. 
While these are all possible square lattices for B < 32, 
this raises the question of whether there are other sym- 
metric minimizers and what are their symmetries, if any. 

It seems that using various spliced initial configura- 
tions, the square symmetry can be relaxed to produce 
other rectangular shapes. We have already seen two of 
these at B = 5 and B = 6, which are not global minima 
and at B = 7 and B = 9, which are the global minima 
for these charges. For all four cases, the short edge was 
the shortest possible, but for B = 10, 12, and 14 a more 
squarelike rectangular lattice of 4 x 5, 4x6, and 4x7 half 
lumps, respectively, exists. These are all global minima; 
they were all obtained using spliced initial configurations. 

The global minima for B = 11 and B = 13 have the 
shape of a distorted rectangle, where the distortion con- 
sists of an extra two half lumps in one corner. These were 
both achieved by using N ~ 801, h = 0.05 initial states 
(as always, these were later coarsened and embedded into 
the standard lattice). Of these, we show the minimizcr 
B = 13 in Fig. \ ^gj\ as an example. 

As the charge increases, finding the global minimum - 
or what can reasonably be believed to be one - becomes 
increasingly difficult. It is likely that the (possibly dis- 
torted) rectangular shapes prevail as global minima for 
charges above B — 14, but we have not been able to as- 
certain this. As an example of how difficult it is to find 
the global minimum, consider charge B = 9. There is a 
local minimum of the distorted rectangle type with two 
extra half lumps in one corner. To converge into this 
configuration, an initial configuration of charges B = 1 
and B = 8 is needed, but it is not enough to just put 
these side by side. They must be placed at very precise 
locations in the lattice. Although this arrangement is 
natural, place the B = 1 solution outside the B — 8 solu- 
tion along one of its diagonals, it is very difficult to guess 
the correct arrangement for higher charges. Moreover, 
at higher charges, there is an ever increasing number of 
different combinations to choose from. 

In addition to the global and polygon-shaped minima, 
we found some other local minima. Most of these have 
energies between the global minimum and the polygon 
and have various (or no obvious) symmetries, like the 
distorted rectangle of B — 9 mentioned above. The one 
notable exception is the lucky cloverleaf at B = 14, whose 
energy exceeds that of polygon and has an obvious sym- 
metry. It is depicted in Fig. | ^h)| 

We now briefly return to the case B = 6. Since it is 
the last one where the polygon seems to be the global 
minimum, it deserves special attention: how can we be 



certain some other configuration does not give a lower 
energy? Since it seems conceivable that a rectangular 
lattice of 3 x 4 half lumps could exist, we tried several 
methods of producing it: we split the rectangular B = 12 
minimizer in two, keeping just one half; spliced together 
all possible combinations of two lower charge solutions 
side by side and even spliced together two B = 1 and 
two B — 2 solitons in a crosslike configuration. They all 
relax either to the polygon or the 2x6 rectangle. Of 
course, it is possible, that the rectangles reported above 
for the higher charge solitons might not be the true global 
minima, but we believe they are: any other rectangle 
would necessarily have a bigger difference between its 
long and short sides. Therefore, we did not attempt to 
find any other rectangles. 

Finally, we summarize the charges, normalized energies 
per B and shapes of our global minimizers in Table HI We 
also conjecture an exponential decay law for E(B)/(4irB) 
of the global minimizer which is approaching the value 
1.52 when B increases without limit. This behavior, 
along with a least squares fit of the form a + be cB , is 
displayed in FigJS] 



B 


Global minimum E/B 


Polygon E/B 


Shape 


1 


1.707 


1.707 


pair of half lumps 


2 


1.591 


1.591 


square 


3 


1.565 


1.565 


regular polygon 


4 


1.556 


1.556 


regular polygon 


5 


1.551 


1.551 


regular polygon 


6 


1.548 


1.548 


regular polygon 


7 


1.547 


1.547 


rectangle 


8 


1.544 


1.546 


square 


9 


1.543 


1.545 


rectangle 


10 


1.541 


1.545 


rectangle 


11 


1.540 


1.544 


distorted rectangle 


12 


1.537 


1.544 


rectangle 


13 


1.537 


1.544 


distorted rectangle 


14 


1.535 


1.544 


rectangle 


18 


1.532 


1.543 


square 


32 


1.526 


1.543 


square 



Table I: Energies of the baby Skyrmions 



VI. CONCLUDING REMARKS 

We have used numerical minimization to show that 
multi-Skyrmions in the easy plane baby Skyrmc model 
exist and are stable. For each charge above 4, we have 
found at least two minima, sometimes three or even 
four. Of these, we have most likely identified the global 
one for most charges in the range [1,14], 18, and 32. 
These minima exhibit various different types of sym- 
metries and it is conceivable that a method exists for 
determining the shape of the minima without finding 
the minima themselves. However, this remains an open 
question. The prevalence of rectangles or slightly dis- 
torted rectangles in this work might provide a way to 
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-0.04 
-0.28 



■ -U.76 

'-1.00 
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(a) Global minimizer for (b) Global minimizer for 
B = 2. B = 6. 




(c) Global minimizer for (d) Global minimizer for 
B = 7. B = 8. 



• • « • • 
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• • • • • 



•••••••• 
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(e) Global minimizer for (f) Global minimizer for 



answer the question. We also present a conjecture that 
lims^ E(B)/(4irB) w 1.52. 

Like baby Skyrme models with some other potentials, 
we have identified an attractive and repulsive channel for 
two well-separated unit Skyrmions using a dipole approx- 
imation and confirmed its validity by numerical simula- 
tions. 

The rectangular lattices formed by the half lumps are 
reminiscent of condensed matter systems, where (baby) 
Skyrmions have been experimentally observed (see, for 
example, [H|,[l4j]). These are not, however, modeled by 
the baby Skyrme model, so it would be interesting to 
investigate whether the easy plane baby Skyrme model 
could be used as an effective model for some condensed 
matter systems. 



B = 18. 



B = 32. 




••• 



■ 4. 

■ -0.7G 
'-1.00 



• • 7 • . 

• •»«•• 



(g) Global minimizer for (h) Local minimizer for 



B = 13. 



B = 14. 



Figure 4: Sample minimizers for higher degrees. All 
plots are contour plots of 4>3 . 



Figure 5: Normalized energy per AttB versus B of those 
configurations which we believe are global minima 
(crosses), the polygon-shaped solitons (circles) and a 
least squares fit of the form a + be cB (dashed curve). 
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